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We report results of a search for the charmless two-body baryonic decays — > pp, AA, 
and pK based on the analysis of a 140 fb~^ data sample. We set 90% confidence level upper 

limits on their branching fractions; ^ pp) < 4.1 xlO"^, S(-B° AA) < 6.9 xlO"^, and 

B{B+ ^pA) < 4.9 xlO"'^. 



PACS numbers: 13.25.Hw 

Motivated by the recent observation of the two-body 
baryonic decay —>■ A+p 0], we search for charmless 
two-body baryonic decays, pp, B^ — > AA, and 

5+ pA. Previous work on these channels set upper 
limits on the branching fractions at the O(10~^) level 

Unexpected, large, transverse polarization has recently 
been measured in B 4>K* decay Q- Measuring the 
polarization in two-body baryonic decays has been pro- 
posed Q as a means to shed light on its origin. This 
could provide the key to understanding whether there is 
new physics in charmless B — > VV decays. 

Recent observations of three-body B decays containing 
baryons in the final states suggest that production is 
enhanced at low invariant masses of the baryon system. 
Some theorists suggest that baryon production is favored 
by the reduced energy release on the baryon side . The 
rates for three-body baryonic decays are large because 
the baryonic daughters can form more readily when an 
energetic meson is recoiling against the di-baryon system. 
In this paper we study the complementary suppression 
of two-body baryonic decays. The results presented here 
are based on a 140 fb~^ data set, corresponding to 152 
million BB pairs, collected with the Belle detector at 
KEKB 0, an asymmetric e+e" collider operating at the 
T(4S) resonance 0. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a three-layer silicon vertex de- 
tector (SVD), a 50-layer central drift chamber (CDC), 
an array of aerogel threshold Cherenkov counters (ACC), 
time-of- flight scintillation counters (TOF), and an array 
of CsI(Tl) crystals, all located inside a superconducting 
solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return located outside of the coil is instrumented to 



detect mesons and to identify muons (KLM). The 
detector is described in detail elsewhere 0- 

Tracking information is collected by the SVD and the 
CDC. For each primary charged track, the impact pa- 
rameter relative to the run-by-run interaction point (IP) 
is required to be within 2 cm along the z axis (aligned 
opposite the positron beam) and within 0.05 cm in the 
plane transverse to this axis. 

Particle identification (PID) for protons, kaons and 
pions is determined from the CDC specific ionization 
[dE/dx), the pulse-height information from the ACC 
and the timing information from the TOF. Proton can- 
didates are selected based on p/ K/tt likelihood functions 
obtained from the PID system. The selection criteria 
are Lp/{Lp + Lk)> 0.6 and Lp/{Lp + L.,^)> 0.6, where 
Lp, Lki and L^^ represent the proton, kaon, and pion 
likelihoods, respectively. The proton detection efficiency 
is 70%, determined from A pn^ decays in the same 
momentum range as B'^ — > pp and B^ pA. The cor- 
responding misidentification rates for charged kaons and 
pions are 11% and 4%, determined using kaons and pions 
from the decay chain D*^ D'^tt^ {K~^Tr~)Tr~ in the 
above momentum range. 

Candidate A baryons are reconstructed in the p7r~ de- 
cay channel and are selected using the following require- 
ments: 1) the distance between two A daughter tracks 
must be less than 12.9 cm along the z axis and greater 
than 0.008 cm in the plane transverse to this axis; 2) the 
fiight length of the A candidates must be greater than 
0.22 cm in the plane transverse to the z axis; 3) the an- 
gular difference between the A momentum vector and the 
vector to the decay vertex from the IP must be less than 
0.09 rad After application of the above selection cri- 
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teria, the reconstructed invariant mass of the A candidate 
is required to be between 1.111 and 1.121 GeV/c^ (Str). 

Candidate B mesons are reconstructed from the pri- 
mary charged tracks in the pp mode; from one pri- 
mary charged track and one selected A candidate in the 
pA mode; from two selected A candidates in the AA 
mode. In these three modes, the B meson candidates 
are selected using two kinematic variables defined in the 
T{AS) center-of-mass (CM) frame: the beam-energy con- 
strained mass, Mbc = \/£'boam ~ ' ^'^'^ energy dif- 
ference, AE = Eg — iS'beam, whcrc p*g and Eg are the 
momentum and energy of the B candidate and E'boam 
is the beam energy. The signal region is defined as 
5.27GeV/c2 < Mbc < 5.29GeV/c2 and \AE\ < 0.05 
GeV. The sideband region, used for background determi- 
nation, is defined as 5.2GeV/c2 < Mbc < 5.26GeV/c2 
and \AE\ < 0.2 GeV. 

The decays K^tt^ and _B" — > tt+tt^ do not 

contribute measurably to the background. This is due 
to their small branching fractions, the small proton fake 
rates, and the shift of the B candidate out of the AE 
signal region. Background from other B decays, studied 
using Monte Carlo (MC) simulation is found to be 
negligible. The main background comes from continuum 
events, especially e+e^ qq {q ~ u, d, s). 

We use the same flavor tagging algorithm used in time 
dependent CP violation measurements to make contin- 
uum suppression more effective. Charged leptons, pions, 
and kaons that are not associated with the reconstructed 
pp, AA, or pA decay are used to identify the flavor of the 
accompanying B meson. The algorithm in Ref. 11] intro- 
duces that the parameter r is an event-by-evcnt dilution 
factor ranging from r = for no flavor discrimination 
to r = 1 for unambiguous flavor assignment. The sig- 
nal and background samples are separated into a group 
of lower tagging quality, < r < 0.75 (ri region) and 
one of higher tagging quality, 0.75 < r < 1 (r2 region); 
these are treated separately in the following continuum 
background rejection method. 

For continuum suppression we use cos^t, which is the 
cosine of the angle between the direction of the primary 
proton (for pp and pA modes) and the thrust axis of 
the non-candidate tracks and showers. For the AA mode, 
it is the cosine of the angle between the direction of the A 
candidate and the thrust axis of the non-candidate tracks 
and showers. We preselect events by requiring cos 9t to 
be less than 0.9. 

We use seven variables to characterize the event topol- 
ogy: five modified Fox- Wolfram moments ll^ll 
and I cos 9t \ ■ We define a Fisher discriminant [15l , 
which is the sum of these seven variables with coefficients 
optimized to separate signal and background events. 
Probability density functions (PDFs) for the signal and 
for background as functions of and of the cosine of the 
polar angle {6b) of the B candidate's flight direction are 
obtained. An additional variable, Az, is used: the decay 



vertex flts are performed using the charged tracks from 
each of the B decays. The distance between the two B 
decay vertices in the z direction, Az, is obtained. For 
BB events the Az distribution is wider than for random 
combinations of proton jets from the continuum, which 
are studied in the MC and sideband data. We use a dou- 
ble Gaussian function for Az as the fltted PDF if \Az\ < 
0.2 cm. For events outside this Az range, we do not use 
the Az PDF for further calculations. 

The signal and background likelihoods, Cs and Cb, are 
formed from the product of the PDFs for T and cos 6b 
and, for the pp mode only, Az. We demand that TZ — 
^s/{^S+^b) exceeds 0.8 (0.85 for the mode) in the ri 
region and 0.2 in the r2 region. This criterion is chosen to 
optimize the flgure-of-merit {Ns/\/Nb) calculated using 
both MC and sideband data samples, where Ns and Nb 
are predicted signal and background yields, respectively. 
We optimize the TZ requirement by assuming that the 
branching fractions for i?° — > pp, AA, and 5+ — > pA are 
1.0 X 10^^. We obtain the expected signal yields from 
the product of the total number of BB events, the signal 
efficiency from MC and the assumed branching fractions. 
The predicted background yields are obtained by fitting 
the AE distribution for the sideband data. 

We compare the sensitivities with and without the in- 
clusion of the Az PDF in the pp mode. With the Az PDF 
included, the figure-of-mcrit improves by 26% in the ri 
region and by 28% in the r2 region. The TZ distributions 
for events in the two regions of r are shown in Fig. ^ 
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FIG. 1: Distributions of TZ for events in the ri and r2 re- 
gions. The open histogram is £s (signal MC) and the shaded 
histogram is jCb (sideband data). The arrows indicate the 
requirement positions. 

The signal efficiency for each mode, after application 
of all the selection criteria, is determined by MC sim- 
ulation and itemized in Table ^ The systematic error 
in this efficiency arises from tracking efficiency (2.0%- 
4.7%), proton identiflcation (0.7%-0.8% per proton, mea- 
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sured from a study of A decays), A selection (2.5% per 
A, including all of the optimized A selection require- 
ments from a study of A decays) , and the likelihood ra- 
tio requirement in the two flavor-tagged regions (3.0%- 
4.6% in the ri region and 4.4%-7.7% in the r2 region, 
based on a study of D'^n^ {K^'it^)tt^ and 

W D+TT^ (i4:-7r+7r+)7r-). By comparing MC 
and data, corrections are made for the effect of pro- 
ton identification (—4.0%, from a study of A decays), 
the likelihood ratio requirement (—4.5% to 0. 8% in the 
ri region and 1.5% to 3.2% in the r2 region, based on 
a study of B- D°Tr- {K-tt+)tt- and B° 
D~^TT~ — !■ (i4r~7r+7r+)7r^), and polarization (—2.2% for 
the B° AA mode, due to the correlation of the spins 
of (p. A) and (p, A)). 

The total systematic uncertainties are 5.1%, 9.2% and 
5.7% for the pp, AA and pA mode, respectively. 

By fitting the AE sideband data (0.05 GeV < |AE| 
< 0.2 GeV), projected on the Mbc axis (5.2GeV/c2 < 
Mbc < 5.3GeV/c^), we obtain the background shape of 
Mbc. By modeling the Mbc distribution with an AR- 
GUS 16] function and assuming a linear shape for AE, 
we obtain the predicted background by scaling the num- 
ber of data events in the sideband region to the signal 
region. 

Since few signal candidates are found (Fig. [21 , we de- 
termine the 90% confidence level upper limits for the 
branching fractions by an extension of the Feldman- 
Cousins method The final results are listed in Ta- 
ble D 
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FIG. 2: Distributions of Mbc and A_B_for (a)(d) _B° -> pp, 
(b)(e) B+ pA and (c)(f) B° AA candidates. We re- 
quire the AE signal selection to plot the Mbc distribution 
and vice versa. The dashed lines represent the fitted back- 
ground curves. 



The results reported here use 4.8 times more BB pairs 



TABLE I: Summary of the B° ^ pp, AA and B+ ^ pA 
search, where e is the reconstruction efficiency, Nobs is the 
observed number of events in the signal region, A^^Ip is the 
expected background in the signal region and BF is the 90% 
confidence level upper limit for the branching fractions. 



Mode 




£ [%] Nots 


' exp 


BF [10-^] 




PP 


20.28 ± 1.03 17 


14.9±3.4 


< 4.1 


SO ^ 


AA 


4.32 ± 0.40 2 


1.2±0.3 


< 6.9 


B+ - 


pA 


9.22 ± 0.53 5 


3.5±1.1 


< 4.9 



than the previous analysis j^- Upper limits on the 
branching fractions at 90% confidence level are 4.1 x 
10"^ 6.9 X 10~^ and 4.9 x 10""^ for the pp, AA and pA 
mode, respectively. The upper limits are reduced by fac- 
tors of 2.9, 1.9 and 4.5. These results are consistent with 
constraints obtained by the other collaborations 

[11, 

and 

indicate that B decays to charmless di-baryon systems 
are suppressed. 

In conclusion, we find no significant signal for the 
modes 5° pp, AA and i?+ — > pA in 152 miUion BB 
events. 
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